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Abstract 

Thermal properties of single species nucleon matter are investigated assuming a simple form of 
the nucleon-nucleon interaction. The nucleons are placed on a cubic lattice, hopping from site 
to site and interacting through a spin-dependent force, as in the extended, attractive Hubbard 
model. A mean field calculation in the Hartree-Fock Bogoliubov approximation suggests that 
the superfluid ground state generated by strong nucleon pairing undergoes a second-order phase 
transition to a normal state as the temperature increases. The calculation is shown to lead to a 
promising description of the thermal properties of low-density neutron matter. A possibility of a 
density wave phase is also examined. 

PACS numbers: 26.60. +c, 21.60.Jz, 21.65. +f 
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I. INTRODUCTION 



Nuclear excitations are complicated dynamical phenomena, depending on the detailed 
structure of the individual nucleus, and must be examined based on the specific structure of 
the nucleus, such as whether it is closed or open shell. As the excitation energy gets higher, 
however, the excitations depend less on specific nuclear structure and start to exhibit more 
common features among (heavy) nuclei. These features are expected to be reasonably well 
represented by excitations of nuclear matter. Furthermore, dynamics of supernovae and 
neutron stars, which have been of much astronomical interest, are expected to be better 
understood through the study of excitations of neutron matter . The gross features of 
;he thermal properties of nucleon matter have been examined by means of statistical models 

and lattice gas models More realistic descriptions of the thermal properties have 
)een provided through applications of various approaches in the nuclear many-body theories 

E3 Q, S 0, 3- Applications of the traditional nuclear many-body theories regarding 
nucleon matter at zero temperature have been extensive and provide the 

most reliable information on the properties of nucleon matter at low temperatures. 

Previously, one of us (R. S.) has collaborated on a Monte Carlo calculation of nuclear 

n 

matter on a lattice [1J] , which provides a new framework for studying the thermal properties 
of nucleon matter. Though the computational space was small and the nucleon-nucleon 
interaction was simple, the calculation has proven to be of much promise, demonstrating the 
occurrence of a phase transition around 15 MeV using the parameters adjusted to reproduce 
the saturation properties. A similar calculation of nucleon matter of a single species has 
also been initiated in the same work, using the same form of the Hamiltonian. Though a 
phase transition appeared to take place at a few MeV, the evidence for it was not quite solid 
owing to statistical fluctuations, which are enhanced at low temperature (a sign problem). 
The phase transitions may correspond to those expected through paired nucleons (Cooper 
pairs) in nucleon matter jll . 

In order to gain a better understanding of the possible, latter phase transition, we apply 
in this work the analytic means of a mean field approach to the problem in the same lattice 
formulation. From this work, we do not expect to be able to draw precise quantitative 
conclusions, but rather we will try to learn the nature of the phase transition at a semi- 
quantitative level. For this purpose, we take the thermodynamical (infinite volume) limit 
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for numerical results, so as to obtain a clear signal of the phase transition. The mean 
field results in this work will also serve as a reference for the more extensive Monte Carlo 
calculation that we are currently carrying out. 

Our Hamiltonian for single-species nucleon matter turns out to be an extended, attractive 
Hubbard model, which has been studied as a simple model of high-temperature supercon- 
ductivity The mean field calculation shows that the low-temperature, low-density state 
is a superfluid state and undergoes a continuous (second-order) phase transition to the nor- 
mal state as the temperature and/or density increases. As the Hamiltonian is not yet fully 
realistic and the values of the interaction parameters are uncertain, our results are not quite 
comparable to those for neutron matter, except perhaps at a very low density. But we 
demonstrate that the approach is promising for the study of low-density neutron matter. 

Furthermore, we find that a density-wave state coexists with the superfluid state, suggest- 
ing that the state of neutron matter may be more complicated than the simple description 
of a superfluid state as it is often characterized. 

The outline of this work is as follows: After the introduction in Sec. I, the Hamiltonian 
and its discretized form in the coordinate space are presented and are identified as an 
extended Hubbard model in Sec. II. The Hartree-Fock Bogoliubov approximation is applied 
and the Hamiltonian is diagonalized in Sec. III. Thermodynamical properties numerically 
calculated are shown in Sec. IV, an attempt to apply our calculation to the problem of low- 
density neutron matter is discussed in Sec. V, and the possibility of a density wave phase is 
examined in Sec. VI. Discussions and conclusion are presented in Sec. VII. 



II. SIMPLIFIED HAMILTONIAN AND EXTENDED HUBBARD MODEL 

The Hamiltonian consists of the kinetic and potential terms K and V, respectively: 
H = K + V 

h 2 



~E/ dr^ T (r)V 2 ^ T (r) 

^ (TT 



+ \ E / dr / dr ' #r(*)#v(*Wr - r0^v'(r')Vw(r) , (1) 

acr' 

tt' 

where is the nucleon mass, and o = ±1/2 and r = ±1/2 are the spin (f or |) and the 
isospin (p or n), respectively. ^ T (r) and Vw( r ) are the creation and annihilation operators 
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of the nucleoli, with the spi n a and isospin r at the position r. As in the previous Monte 



14j . we include only the central and spin-exchange interactions, V c 



Carlo lattice calculation 
and V a , respectively: 

V(t - r') = V c (r - r') + V (r - r')tr ■ a' . (2) 
V c and V a are taken to consist of on-site and next-neighbor interactions, 



V c (r - r') = V®6 (r - r') + V^ [V 2 r 6 (r - r' 
V a (r - r') = Vj 0) 5 (r - r') + V} 2) [V 2 r 5 (r - r' 



(3) 



where V^ terms can be written explicitly exhibiting their hermiticity. 

As physics of the lattice description is more apparent in the coordinate space, we consider 
the discretized coordinate with an inter-nucleon spacing a in the cubic lattice with the torus 
boundary conditions. We thus focus our interest on the physics of the spatial separation 
greater than a in each direction, or of the momentum component roughly between 71/ a 
and —71 /a, by eliminating (or integrating out) the physics of the shorter distance. The 
discretization corresponds to 



where % denotes a lattice site specified by rij with its component ranging [—aN/2, aN/2]. 
Here, N is the number of sites in each spatial direction. Note that the creation and anni- 
hilation operators, c\ ur and c i(TT , have no dimension as defined. We also apply the identity 



i 

where i denotes the spatial components (x,y,z), and the Greek indices denote the compo- 
nents of the Pauli spin matrix, 1 or 2. 

In this work, we study the simplified case of nucleon matter consisting of a single nucleon 
species, such as neutron matter. The Hamiltonian Eq. ((TJ is then expressed in a spatially 
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discretized form: 



{i,j)cF io- i 

(4) 

(mV (ij)o- 

where denotes the pairs of next-neighbor sites, and t is the hopping (kinetic energy) 
parameter defined as 



2mjva 2 

Here, the potential parameters U and V^'s are expressed in terms of linear combinations of 
V^'s and 1/ (2), s. The Hamiltonian Eq. (gj) is now in the form of an extended Hubbard 
model, which is the Hubbard model with the on-site spin-pairing interaction of the U term, 
modified by the next-neighbor interaction of the V terms. As the U value will be taken to 
be negative, our model is an extended, attractive Hubbard model. The repulsive Hubbard 
model has been well studied in condensed matter physics as a model of strongly correlated 
electron systems Q], but the attractive model is generally less studied. In recent years, 
however, the extended, attractive Hubbard model has drawn much attention as the model 
describing the essential features of high-temperature superconductivity [15(. Note that the 
extended, attractive (negative-?/) Hubbard model used in condensed matter physics, how- 
ever, is usually of the simpler form shown below and has no 6t term as a part of the kinetic 
energy \v\ . 

When the spin- dependent next-neighbor interaction is taken to be small and negligible, 

K (2) = , (6) 

the Hamiltonian is simplified: 

H = ~t ^iaCja + 6t ^ 6 tcia + U ^ C^C;^ + V ^ c\ a C ia c) a ,C j(jl , (7) 

(i,j)a icr i (i,j)cra' 



where 



V = • (8) 
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Apart from the lattice spacing a, the Hamiltonian of Eq. (JJJ) now describes dynamics with 
two parameters. The Hamiltonian Eq. (j3J) [and thus also Eq. (JJJ)] possesses an underlying 
particle-hole symmetry, which affects thermodynamical properties as discussed in Sec. VI. 
The symmetry is not an explicit property in our original Hamiltonian, Eqs. (l)-(3). We 
elaborate on the symmetry in Appendix A. 



III. HARTREE-FOCK BOGOLIUBOV APPROXIMATION AND GAP EQUA- 
TIONS 

We now apply the mean field method in the Hartree-Fock Bogoliubov approximation 
UM. Here, we expect effects of the U term to dominate the thermal properties of the single 

nnn 

species matter as in the standard BCS description |19|, |2C| |21[ , but we also wish to treat 



their single-particle aspects in the Hartree-Fock approximation on the same footing. Since 
the method is well known, we limit the description of the formalism to the key steps that 
are specifically relevant to our calculation. 

The nature of the mean field approximation is apparent in the spatial representation. 
Through the application of the Wick theorem, the decoupling scheme for the U term is 

where 

A* = (c^Cii) = -{ciiCi}) , 

n i] = (c! T c iT ) and n i{ = (cj^) . 

Here, ( • • • ) denotes the expectation value in the BCS-like ground state, which is to be 
determined in a self-consistent way. The order parameters Aj's are related to the local 
density of the condensate of nucleon pairs, while n^'s are the average number of nucleons. 
Note that the inclusion of rii as a variational parameter distinguishes the present treatment 
from the standard BCS Q|. For the V terms, we have 

t A At . _ At A. . , „. At 



C icr C iv C 



In the following, Aj and n^'s will be assumed to be independent of the site i, or global, 
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and will be formally treated as the variational parameters: 

A ~ 2UAi = 2UA* , 
n ~ 2n iT = 2riii . 

Our A is denned to be the gap energy itself, with the dimension of energy in the unit of 2U, 
and carries an extra factor of 2 in comparison to the often-used A. We also note that we 
ignore the V-term contribution to the hopping term, as they merely change somewhat the 
strength of the hopping term and of the constant part of the energy, without affecting the 
physics of the phase transition. The mean field approximation then yields in a cubic lattice 
with six neighboring sites 



where U = U + 24V. 

For convenience, we will carry out the mean field calculations in the momentum 
space. The momentum representation is introduced by discretizing the momentum as 
p — > 27rk/(iVa), with each component of k being an integer, ranging [— N/2, N/2]. Note 
that we now have r • p — > 27rnj • k/N. The coordinate and momentum representations of the 
operator c are related through the Fourier transformations 



and similarly for the nucleon creation operator, c' . Note that the discretized orthonormality 
relation is 







(10) 
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With the chemical potential //, H — jj,N 3 n is in the momentum space 




(11) 
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where 

ji = ft - (Qt + nU/2) , (12) 

and 

e k = -t exp (2mk • e/N) = -2t ^ cos ( 27rk j/ N ) ( 13 ) 

e j=x,y,z 

is a part of the kinetic energy of a quasi-particle expressed in terms of the unit vector e 
showing a next-neighbor site. Note that e_ k = e k and X^k e k = 0- We see that in the 
Hamiltonian, V appears only as U = U + 24V, and merely shifts the chemical potential and 
the total energy: It does not actively participate in the generation of the phase transition. 
As is well known, the spin density s, 

s = ^ = ^Jp _ 6 lA> = Z^At " a -ki a -ki> , (14) 

i k 

is conserved in the mean field approach, while n, 

n= (h) = — ^(44x) = ^3 J2(cLc ka ) (15) 

iff kcr 

(and thus the total nucleon number), is not, as the spin number operator s commutes with 
H of Eq. (0), but the number density operator n does not. We remedy this problem by 
the standard method of introducing a Lagrange multiplier corresponding to the chemical 
potential /x, by subtracting a term fiN 3 n = xx SoA'ct fr° m the if. /x will be adjusted to 
achieve the desired value of the conjugate parameter, n. The formalism is thus essentially 
the canonical ensemble method. 

We now diagonalize H — fiN 3 n by the use of the Bogoliubov-Valatin transformation, 

/3l k _ = w k cl ki + v k c kT , 

where u k and t> k are taken to be real and are given by 

•i 1 f, , ek~P>\ ,2 1 A e k -/x 
u \c — — I 1 H ^ — I and v k = - ( 1 



2 V i? k / 2 y i? k 

satisfying M k + f k = 1. H — fiN 3 h is now expressed as that of a system of free quasi-particles: 

H - fiN 3 n = £kALAA + N 3 (E GS - fin) , (17) 

kA=± 
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with the energy of a quasi-particle 



E k = y/(e k - ft) 2 + A 2 /4 (18) 
and the ground-state energy of the system 

E GS = ~\(^ + n 2 u\ -A- 3 Y, E *-» + » n - ( 19 ) 

' k 

Equation shows that /3 kA 's obey anti-commutation relations and that the quasi-particles 
are fermions. Furthermore, Eq. (|17j) implies that they form a system of free fermions. As a 
consequence of the thermal average, the internal energy is then given by 

E = (H) = J2 ^kA + N 3 E GS , (20) 

kA 

with the energy per lattice site being E/N 3 . Here, n kA is the momentum distribution of the 
quasi-particles, n kA = (/3,J, A /3 kA ), and is determined by the requirement that the free energy 
F introduced below is minimized by a variation of n kA , SF/dn^x = for A = ±. We obtain 

n k = n k+ = n k _ = [exp (E k /T) + l]" 1 , (21) 

which have the limiting values, n k —>■ and — >• 1/2 as T — >• + and — > oo, respectively. 
Note that throughout this work, we denote the temperature T in the unit of the Boltzmann 
constant. 

By combining with the entropy S, the (Helmholtz) free energy is expressed as 

F(T, a; A, n) = E-TS 

= ^£ k ^kA + iV 3 £ GS 

kA 

+T^[n kA lnn kA + (l-n kA )ln(l-n kA )] , (22) 

kA 

depending on a through the k sum because the spatial volume (aN) 3 depends on a with 
iV fixed. F is a function of A and n. A and /i are determined so as to minimize F for 
variations of A and n, while T and a are fixed. The conditions 



provide the gap equations 

e k - p, 




(2n k - 1) , 

(23) 



from which \x and A are determined. 

IV. THERMODYNAMIC AL PROPERTIES 

We now apply the formalism so far described, to compute various thermodynamical quan- 
tities. For clarifying our presentation, we place some expressions of the thermodynamical 
variables in Appendix B. All numerical results are calculated in the thermodynamical limit 
N — > oo. In the limit, the summation over the discretized momentum space of each com- 
ponent of k ranging in [— N/2, N/2] is replaced by the integral over the first Brillioun zone 
with each component of the momentum p ranging [— ir/a, ir/a]: 

3 pir/a p7r/a pn/a 



jvsXXs) / / / ' i' 



■k I a J — n/a J — ir/a 



A. Potential parameters 

We apply the parameter values used in the previous Monte Carlo lattice calculation for 



nuclear matter 



V c {0) = -181.5 MeV fm 3 



V c {2) = 37.8 MeV fm 5 , 

(0) 3 ( 24 ) 

Vr = -31.25 MeV fm 3 , 

K (2) = , 

with the lattice spacing a = 1.842 fm. These parameter values give 

t = 6.11 MeV , 
U = -24.74 MeV , 
V = 0.89 MeV . 

The parameter values of Eq. ()24|) were chosen in Ref. Q| so as to reproduce the saturation 
density and energy of nuclear matter on a finite 4x4x4 lattice for the same Hamiltonian as 

10 



ours, Eqs. (HJ-©- Our Hamiltonian has no explicit r-dependent term, and the parameter 
values effectively include the strong neutron-proton interactions for nuclear matter. The use 
of the parameters is thus not quite adequate as a realistic description of the nucleon matter 
of single species, such as the neutron matter. Furthermore, finite lattice volume effects make 
the thermodynamical limit (N — > oo) calculation different from the finite volume calculation. 
For comparison purposes with the previous and future Monte Carlo calculations, however, 
we use the above parameter values except when the U dependence of A is examined. Any 
conclusion that we could draw from the numerical results in this section is then qualitative. 



B. Gap parameter A 

Equations ()23j) determine A and /i. Figure ^ illustrates A as a function of the temperature 
T for n = 0.5 (one-quarter filling), 1.0 (one-half filling), and 1.5 (three-quarter filling). In 
the figure, we see that A vanishes at T = T c , 

0.66* or 4.0 MeV for n = 1.0 , 
0.55t or 3.3 MeV for n = 0.5 and 1.5 . 
The temperature dependence of A is the same for n = 0.5 and 1.5. This is a consequence of 
the symmetry with respect to n — 1 and is discussed further in Sec. IV D and Appendix A. 
Figure ^ also shows 

A(T = 0) ~ (2 - 2.5)t ~ 12 — 15 MeV . 

The explicit values of T c and A(T = 0) depend sensitively on the parameter values as 
discussed below, but T c and A(T = 0) satisfy 

A(T = 0) ~ 3.6T C , 

which is the well-known relation at the weak-coupling limit except for the latter to have 
a slightly smaller coefficient 3.54 [21]. In comparison to the weak limit, our calculation 
thus somewhat underestimates T c in relation to A(T = 0). Note that different mean field 
calculations have been reported to yield the coefficient smaller than the weak-limit value 
jiol ] and also even much larger jjj than ours. 
Near T c , the gap equations Eqs. (|2*3*jl yield 

(T - T c f for T < T c , 



A oc 



for T > T c , 
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with (3 ~ 0.45. Note that the well-known mean field value of the critical exponent (3 in the 
simple BCS theory is 1/2 3]- 

The physics of the phase transition depends on the strength of the potential parameter 
U. Figure 121 illustrates how sensitively the value of A depends on U. We see that A — ► 
as U — > 0. As is well known, A does not vanish for a finite U. In fact, Eq. (|23|) yields the 
well-known dependence of A on \U\ for U — > 0: 

A -> Ae-^/l^l , 

where A and -B are constant and independent of U. 



C. Second-order (continuous) phase transition 



The temperature dependence of A in Fig. ^ is a well-known dependence of the order 
parameter for a second-order phase transition. The variation of A as shown in Fig. ^implies 
that the phase transition takes place from a superfluid state generated by spin pairing to 
the normal state, as the temperature increases. The features of the second-order phase 
transition are clearly seen in the temperature dependence of the thermodynamic quantities 
expressed in terms of the temperature derivatives of the free energy F in a successive order. 
We consider the internal energy E, the entropy S, and the heat capacity C v , 

~d{F/T)~ 



E = 
S = 

c v 



dT 



OF 



dT 



(25) 



-T 



OS 
df 



-T 



d 2 F 
-Qfl 



The temperature dependence of the quantities is calculated using Eqs. (|2~Uj) . (|2*2*|) . ()B1|) . and 
flB2|) . and is illustrated in Figs. 01 and |U E and S are continuous at the critical temperature 
T c as seen in Fig. 01 while C v has a jump at T c as in Fig. 0] These behaviors demonstrate 
the generic features of the second-order phase transition. The amount of the discontinuity 
in C v at T c , AC V , is relative to C v of the normal phase, 



AC„/ C v (normal) 



1.74 for n = 1.0, 

1.43 for n = 0.5 and 1.5 
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while the BCS mean field value is 12/7^(3) ~ 1.43, independent of n [21]. (( is the zeta 
function.) 

The thermal quantities involving volume derivatives form a set of quantities similar to the 
temperature derivatives. We consider the pressure P and the isothermal compressibility k^. 
Here, following the common practice in nuclear physics, we examine the incompressibility 
K = 9/ (^tp), defined in terms of the density p = V / (aN) 3 = n/a 3 with the spatial volume 
V = (aN) 3 . A volume derivative is then a derivative with respect to the lattice spacing a. A 
derivative with respect to a requires, however, the knowledge of a dependence of U and V, 
that is, their renormalization flow when a is varied. In this work, for simplicity, we assume 
that their a dependence is small, at least around the value of a we use. P and K are written 
as 

~dF~ 

■> 

T,a 



P = - 



dV 



V 

K = -9- 
P 



~dP~ 




-9* 


~d 2 F~ 




_dV_ 


T,a 


P 


dv 2 _ 


T,a 



(26) 
(27) 



P is calculated using Eq. (jB4J) and K is obtained numerically from the temperature de- 
pendence of P. The temperature dependence of P and K confirms the second-order phase 
transition, as seen in Fig. |3J 

There are many other quantities that describe thermodynamic properties of the system 
described by our Hamiltonian, but they are either related to the quantities already shown, 
or their features depend strongly on the explicit form of the Hamiltonian. We thus do not 
show them in this exploratory work. For example, the temperature dependence of double 
occupancy per site, 

D ^ ^E^MA) > (28) 

i 

provides the amount of the spin pairing that participates in the phase transition. But D, 
as well as the kinetic energy per site, is related to A and n in the mean field approximation 
23| and provides no new information, as shown in Appendices A and B. 



D. Particle-hole symmetry and phase diagram 



For illustrative purposes, however, we show the density dependence of /x, D, E, and KE 
at T = 0, in Fig. EJ As noted above, Fig. [T] suggests a symmetric dependence of A on 
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the density n, which is more clearly illustrated in Fig. |3 The symmetry is generated as a 
consequence of the particle-hole symmetry, as seen from the fact that Eq. is invariant 
under the particle-hole conjugation though the symmetry is implicit in our Hamiltonian, 
Eqs. (jlj and (J7j). There is a group of Hamiltonians, in which the particle-hole symmetry 
is implicit, yet yielding (in the mean field results) the energy spectrum of the system with 
the explicit symmetry. The Hamiltonians do yield different behaviors of thermodynamic 
variables. Figure El is an example. We elaborate on the issue of this symmetry in Appendix 
A. 

Some thermodynamic variables are made to exhibit this symmetry explicitly by modifying 
the Hamiltonian to possess the explicit symmetry. Whether it is explicit or not, however, 
the thermodynamic properties obtained from the Hamiltonian are affected by the symmetry. 
The symmetry is a consequence of our computational method using the lattice configuration, 
and it is an artifact. In order to extract physically realistic results, we should therefore stay 
away from the region of the symmetry and should confine ourselves to a small value of n 
by appropriately adjusting the value of the lattice spacing a, so as to simulate the desired 
density of the nucleon matter. We discuss this point again in the following section, where 
we attempt to apply our calculation to a case of low-density neutron matter. 

Combining the variations of the thermodynamic quantities, some of which have been 
presented so far, we obtain the phase diagrams of the present system described in the mean 
field theory. Figure [3 shows the phase diagram in the region of small densities where we 
expect the above-mentioned symmetry to be generating less distortion. 

V. LOW-DENSITY NEUTRON MATTER 

In the previous section, we have used the parameter values most appropriate as a descrip- 
tion of nuclear matter and have examined the nature of the single species nucleon matter 
described by our model. In this section, we discuss whether our model could be made a 
realistic description of neutron matter. 

First, we have the question of whether our lattice would meet the basic momentum 
requirement imposed by the lattice spacing. A lattice description can be made realistic 
when the lattice spacing is less than the momentum scale of the system. By taking the 
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Fermi momentum of the neutron matter as an estimate of the momentum scale, we have 

n/a>p F . (29) 

Our lattice spacing, a = 1.842 fm, yields the density of p = 1/a 3 ~ 0.160 fm -3 for n = 1.0 
(the lattice space being half full). The Fermi momentum corresponding to this density is 

Pf — 1-68 fm -1 

using p = p F 3 /3iT 2 . The value of p F is practically the same as 

n/a ~ 1.71 fm -1 , 

and thus the lattice with the above lattice spacing is applicable to a density much smaller 
than p ~ 0.160 fm" 3 . Note that the preceding discussion yields that the condition 

7r/3 > n 

meets Eq. fj2~9"j) independently of a. 

Second, there is the question of whether our Hamiltonian is appropriate for a realistic 
description of low-density neutron matter. The nucleon-nucleon 1 Sq phase shift is much 
greater than the nucleon-nucleon phase shifts of other states below the laboratory energy 
Ei a b ~ 100 MeV. The nucleon momentum in the center-of-mass coordinate system p cm 
corresponding to this Ei ab is about 1.2 fm _1 through E iab = \{jp 2 cm l'lrn N ), and is smaller 
than the above ir/a = Pcutoff — 1-7 fm" 1 . We thus infer that our Hamiltonian form of the 
S- wave should be reasonable for neutron matter of a density less than 0.17 fm -3 , which 
corresponds to the Fermi momentum of 1.7 fm" 1 . 

As to the parameter values in the Hamiltonian, it would be best to determine them 
for our lattice size from experimental 1 So phase shifts by applying the method of effective 
field theory 0,13- Instead, as an exploratory study, we simply adjust the [/-parameter 
value so as to see whether our approach could come close to other mean field calculations 
of low-density neutron matter in the literature. Figure |H1 illustrates that we could obtain a 
somewhat reasonable density dependence of A by increasing the magnitude of U. We leave 
a more serious determination of the parameters for our future work. 
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VI. DENSITY WAVE PHASE 



Our analysis has so far been strictly based on the mean field approximation applied to 
a spin-pairing phase at the same site. In the approximation, the V term merely shifts the 
effective chemical potential and is inactive in generating the phase transition. As the V 
term represents the pairing of the spin densities at the adjacent sites, such a role may be 
a reasonable one in this phase transistion. Would the V term ever play an active role in 
generating a different phase transition? In this section we briefly examine this possibility. 

The most likely phase in which the V term would play the major role would be a density 
wave phase generated by a coupling of the densities of the opposite spins at the adjacent 
sites. We examine how the V term could generate such a phase transition, again in the mean 
field approximation, and see whether the phase transition would occur with our parameter 
values. 

We follow the Hartree-Fock Bogoliubov approximation with the same decoupling scheme 
as before. For simplicity, however, we ignore the superfluid phase. The density wave phase 
is introduced by making a replacement, 

+ (5/2) cos(27rni • q/N) 
~ n/2 + (5/2) cos(27r ni • q/N) . (30) 

Here, the inhomogeneous order parameter for the density wave depends on the amplitude 5 
and the wave number vector q. In the following, we examine the wave modes in which the 
adjacent sites are the maximum and minimum of the amplitude in a cubic lattice. That is, 
at least one component of 2nq/N is ±n. The number of the non-zero components, or the 
dimension of the density wave d, provides a convenient parameter rj, 

7] = 3 - 2d . (31) 

As in the previous spin-pairing case, n and 5 are treated as independent parameters, and /i, 
the parameter conjugate to n, is adjusted to the desired value of n. 

The steps to the gap equations are similar to the case of the superfluid phase discussed 
in Sec. Ill, and are shown in Appendix C. By minimizing F — fiN 3 n, we obtain 
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by varying n and 5, respectively. Equation (|32|) determines 6 and /i. Furthermore, as 
«k+ < "ftk- and 

C/ + 8r]V < (33) 

for our parameter values, Eq. (|32|) shows that the density waves of all dimensions are expected 
to occur owing to the strong, attractive U. 



VII. DISCUSSION AND CONCLUSION 

The form of the nucleon-nucleon interaction, Eqs. (l)-(3), shows that the nucleon-nucleon 
interaction used in this work is of the S states. As is well known, the realistic nucleon-nucleon 
interaction is highly state-dependent. The relevant nucleon energy of interest to us here is a 
few hundred MeV, corresponding to the Fermi energy region of the nuclear matter density. In 
this energy region, the attractive neutron-neutron interaction is known to be dominated by 
the 3 P2 interaction driven by the spin-orbit force, coupled with the 3 F2 interaction associated 
with the tensor force 27|,|28|. Our Hamiltonian accommodates none of these features of the 
interaction. Our objective as noted in Sec. I is to understand the essential physics associated 
with the thermal properties of nucleon matter, but our finding in this work is limited in this 
sense and is perhaps most applicable to low-density neutron matter. 

There is a serious question of how good the mean field calculation is in our case. A 
mean field approximation ignores most features of particle correlations. As the particle 
correlations are the vital ingredient of critical phenomena, the mean field approximation 
is generally believed to be only of qualitative use, and some times not even qualitative, as 
fluctuations could alter the nature of the phase transition. The situation, however, depends 



on the nature of the problem 22], as the prominent success of the Ginzburg-Landau/BCS 
theory shows j^J, especially at temperatures not too close to the critical one. Our problem 
is in three dimensions, close to the usual upper critical dimension of four under the Ginsburg 



criterion 



221 ] . We are hoping that our calculation, being similar to the BCS theory, is not far 



off, but this remains to be seen. This issue is under further investigation b y in corporating a 



renormalization approach, as has recently been done at zero temperature 30 1. 

Our Hamiltonian has a form similar to that of the Skyrme interaction 3l( (though ours 
is a truncated form.) The Skyrme interaction is one of the effective interactions that are 
phenomenologically introduced to achieve quantitative agreement with experiments, usually 
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by the use of a mean field approximation such as the Hartree-Fock calculations. Though 
it is still not quite realistic, our Hamiltonian has a justification in this sense. A lattice 
calculation such as the previous Monte Carlo lattice calculation [l4j . however, accounts 
for all the complexity of the many-body interaction with no approximation other than the 
numerical, in the lattice framework. The nucleon-nucleon interaction used in it should not 
be then an effective interaction like the Skyrme interaction, but an interaction in free space. 
The parameters U and V are expected to be determined from scattering data through the use 



of effective field theory by extending Liischer's formula [32J, |33( for the large 1 S'o scattering 
length. This issue is presently under investigation [2^. Note that the values of U and V 
will depend on the lattice spacing a. 

In conclusion, our mean field calculation on a lattice with a simple nucleon-nucleon inter- 
action suggests a second-order phase transition taking place at a low temperature in single 
species nucleon matter described by a simple Hamiltonian. Thermodynamic variables show 
their dependence on the temperature and density variations as expected under the phase 
transition in the mean field approach. The transition changes the phase of the matter from 
a superfluid state due to nucleon pairing, to a normal state as the temperature increases and 
the density decrease. This dependence is in qualitative agreement with the findings that 
have been reported in the literature 
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APPENDIX A: THE PARTICLE-HOLE SYMMETRY 



The Hamiltonian of Eq. (J7J) does not explicitly exhibit particle-hole symmetry, but it can 
be modified to do so by adding a single-particle Hamiltonian, 

AH = N 3 U/A- (6t + U/2) ^n ia , (Al) 

where n ia = c\ a c i(T . The new Hamiltonian, 
H' = H + AH 



n jV , (A2) 



(i,j)aa' 

is symmetric under the particle-hole conjugation (or with respect to half filling, n = 1): 



st 

Cj CT C i(T , 



or nj CT <-> 1 — which is n <-> 2 — n with n = ^2 ia .n ia . Under the conjugation, the 
Hamiltonian becomes that of holes, with t — > — t. 

When we repeat the BCS formulation in Sec. Ill, after the Bogoliubov transformation, 
we obtain the diagonalized form of the new Hamiltonian, Eq. (|A2|) . as 

H' = J2K<x + N 3 E' GS (A3) 

kA 

with the quasi-particle energy 



E' k = v /(e k - / x') 2 + A 2 /4 (A4) 

and 

-4 k = /t k+ 
Here, the ground-state energy is 



nl = nL_ = n' k _ = [exp {E'JT) + 1] 1 . (A5) 



1 n 2 - l-\ 1 . 



(A6) 



19 



with 

p! = n-{n- l)U/2 . (A7) 

Equations (fT2j) and ()A7|) show that the \i value is shifted, and Eqs. (fT9j) and (jA6J) tell us 

that the expression of the ground-state energy is altered. The kinetic energy is also changed: 

in comparison to Eq. ()B7|) . We also see that _Eq S and KE' for the new Hamiltonian explicitly 
exhibit the particle-hole symmetry about half filling (n — 1). 
The double occupancy per site for the new Hamiltonian D', 

D> = W* 5>lA^i> = ^ + x ' (A9) 



differs from that for the original Hamiltonian, Eq. ()B6j) . Note that neither D nor D' is 
symmetric about n = 1. 

The gap equations for the new Hamiltonian become 



n - 1 = ^E^(H-i) 



(AlO) 



These gap equations are the same as Eqs. except // is replaced by //'. As we solve the 
gap equations for A and \x (or Jx) for a fixed n, the gap equations for the two Hamiltonians 
yield the same set of A and p,, thus the same T c . Both Hamiltonians thus provide the same 
excitation energy spectrum. 



APPENDIX B: THERMODYNAMIC VARIABLES 

We list here the expressions of the thermodynamic variables, which are used in Sec. IV. 
The entropy is given in terms of n^x in Eq. (}2*Tj) by 

S = t nkA ln ^kA + (1 - n kX ) In (1 - n kA )] . (Bl) 

kA 
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By the use of Eqs. (fTHj). (pHJj) . and (J22I), the heat capacity C v is expressed as 
v d 2 F - 



Cv 



-T 



Q T 2 



1 AA T , , 2 „ / 

u * i, \ 



AAr _ #§\ Sk/r _ AAr 
4 T 2 / 217 



Here, A T is 



At = 



dA 

df ' 



and is numerically calculated from the solution of the gap equations, Eq. (J23|) . 
The pressure P is written in terms of the space volume V = (aN) 3 as 
~dF~ 

T,n 



P 



1 dE\, , n sn \ ( dU 



da 3 



where 



dE k 

da 3 
de k 
da 3 

dU 



dtv n ( dU 



+ - 



da 3 2 \ da 3 



+ 24 



dV 

da 3 



3a 2 V ^ J 



i=x,y,z 



da 3 
dV 



V (o) 



V, 



(2)' 



3 V® 



(2)~ 



9a 3 6a 8 c 
The double occupancy per site D is expressed as 



D 



/V3 2^t Ci T C i! C ^/ - 



AT3 

A 2 n 2 n 
+ -r ~ 



4XP 4 N 3 



The kinetic energy per site KE is written as 
KE e 



iV 3 ^ 



(e k - //) (e k + 6t) 



[2n k (l - -^ k e^ 



+ 6t . 
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APPENDIX C: DENSITY WAVE PHASE 

The Hamiltonian in the coordinate space is reduced to 

h = -t J2 4^ + ( Qt + n u/2) 4 

+(U/2 + A V V)S cos(27r ni • q/N)^ 

ia 

- [{U/A) (n 2 + 5 2 /2) + V(6n 2 + V 5 2 )] N 3 , (CI) 

where the parameter rj depends on the dimension of the density wave d, 77 = 3 — 2d. Here, 
we have used 

(i,j)<T (ij)cr iecr 

~ 6n^cLQ (T + 2^Xcos(27rn i -q/^cLc^ (C2) 

ia ia 

and 

^(n iT + ni;)(n iT + n i4 ) -> ^(n; T + n i; )(n iT + n^) + 25 ^ cos(27rn* • q//V)(n iT + n n ) 
{id) (hi) foV 

+5 2 X cos(27rrii • q/iV) cos(27r(rii + e) ■ q/N) 

ie 

~ (6n 2 + r^ 2 )iV 3 , (C3) 

where the nucleon densities n's are defined in the same way as in Sec. III. In the last step 
of Eq. (|C3|) . we have used the identities 

X cos(27rni • q/N) = 

iV 3 ( C4 ) 
Xcos 2 (27rn,-q/iV) = — 

i 

by applying the discretized orthonormality relation Eq. (jlOj) . 
Using 

S COS ( ^11, • q c\ a C la = ( £ Lq/2 <A+q/2 <x + C^+q/2 <A~q/2 a) , (C5) 

ia ^ ktr 
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)l 









we obtain the momentum space representation 

kcr kcr 

( £t - £ +- + - \ f5N * + ^ 

kcr 

-LsSN 3 + EoN 3 
16 

where the subscripts + and — denote k+q/2 and k— q/2, respectively; that is, c± a 
e± is defined as 

e± = -t^ exp (i2vr(k ± q/2) • e/N) = -2t ^ cos ( 2n ( k i ± <lj/ 2 )/ N ) 

e j=x,y,z 

We have also defined 

5 = 5 (2U + lQr]V) . 
H — fiN 3 h is diagonalized through the transformation 

«k+ = «kC+ CT + ^kC- CT , 
«k- = WkC.^ - v k c +a , 
where «k and t>k are taken to be real and are given by 

satisfying + — 1. Here, 

H — fiN 3 h is now expressed as that of a system of free quasi-particles: 
H - flN 3 h =lYl ^kA«L«kA + N 3 (E GS - H , 

where A = ± and 



2 

kA 



£k± = i[(e + + e--2/i)±6c] 
is the energy of the quasi-particles, and 

E GS = -66/16 - n 2 U/A + un 

is the ground-state energy of the system. 
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FIG. 1: The order parameter A as a function of temperature T in the unit of hopping parameter 
t for the density n = 0.5 and 1.5 (long-dashed curve) and n = 1.0 (solid curve). 



5 2 



— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 


T 1 1 1 1 1 1 1 1 1 1 1 1 


^ \ 
\ \ 

" N \ 

N \ 

\ \ 
X \ 

N \ 
X \ 
\ \ 


n = 1.0 : 


X \ 

! n = 0.5, 1.5 
, , , , i , . . . i . 





-5 -4 -3 -2 -1 

U/t 

FIG. 2: The order parameter A as a function of the on-site potential strength U in the unit of 
hopping parameter t at zero temperature (T = 0) for the density n = 0.5,1.5 (long-dash curve), 
n = 1.0 (solid curve). 
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FIG. 3: The Helmholtz free energy F (solid curve), the internal energy E (long-dashed curve), and 
the entropy of the system S (short-dashed curve) as a function of temperature T in the unit of 
hopping parameter t for the density n = 1.0. The left axis refers to F and E, and the right one is 
for S. 
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FIG. 4: The heat capacity C v for the density n = 0.5 and 1.5 (long-dashed curve) and n = 1.0 
(solid curve), the pressure P for the density n = 1.0 (short-dashed curve), and the incompressibility 
K for the density n = 1.0 (dotted curve) as a function of temperature T in the unit of hopping 
parameter t . 
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FIG. 5: The chemical potential fi (solid curve), the double occupancy D (long-dashed curve), the 
internal energy E (short-dashed curve), and the kinetic energy KE (dotted curve) in the unit of 
hopping parameter t as a function of the density n at zero temperature (T = 0). The left axis 
refers to fi, E, and KE, and the right one is for D. 
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FIG. 6: A in the unit of hopping parameter i as a function of the density n at zero temperature 
(T = 0). 
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FIG. 7: T-p phase diagram. 
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FIG. 8: The gap energy A as a function of the Fermi momentum kF for U = —45.8 MeV (solid 
curve) in comparison to a simple BCS calculation (long-dashed curve) and a more elaborate one 
[l(il | (short-dashed curve). The simple calculation is taken from Ref. j^] [using Eq. (8) with the 
neutron mass of 940 MeV and the neutron-neutron scattering length of -18.8 fm]. A in this figure 



is half the A defined in the text, so that its definition agrees with the one used in 
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